Exponential functionals of Brownian motion and class one Whittaker 

functions 
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Abstract. We consider exponential functionals of a multi-dimensional Brownian motion with 
drift, defined via a collection of linear functionals. We give a characterisation of the Laplace 
transform of their joint law as the unique bounded solution, up to a constant factor, to a 
Schrodinger-type partial differential equation. We derive a similar equation for the probability 
density. We then characterise all diffusions which can be interpreted as having the law of the 
Brownian motion with drift conditioned on the law of its exponential functionals. In the case 
where the family of linear functionals is a set of simple roots, the Laplace transform of the 
joint law of the corresponding exponential functionals can be expressed in terms of a class one 
Whittaker function associated with the corresponding root system. In this setting, we establish 
some basic properties of the corresponding diffusions, which we call Whittaker processes. 



1. Introduction 

Matsumoto and Yor |18j proved that, if (B^(t),t > 0) is a standard one-dimensional Brow- 
nian motion with drift [i, then 



(b^ (t)+ log (j\- 2BM ^ds^ ,i>0 



is a diffusion with generator given by 

where is the Macdonald function. One of the motivations for this work is to identify a class 
of diffusions which should play a similar role in a multi-dimensional generalisation of the theorem 
of Matsumoto and Yor, analagous to the multi-dimensional generalisations of Pitman's '2M — X' 
theorem obtained in [5j 1231 14], 

Note that — K-^. If [i > 0, the diffusion with generator (JT|) can be interpreted as the 

Brownian motion conditioned on the exponential functional A x — e^ 23 ^ ds having a 
certain distribution (a Generalised Inverse Gaussian law). This interpretation is made precise in 
the paper [3]. 

In this paper we consider a Brownian motion flM in W 1 with drift [i, and a collection of linear 
functionals ai, . . . ,ay such that the exponential functionals 

/>oo 

Al = / e-^^ds (i=l,...,d) 
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are almost surely finite. Our aim is to understand which diffusions can arise when we condition 
on the law of 4 M = (A^, . . . , A^). The first step is to understand the law of A^. We show that 
the Laplace transform of A^ satisfies a certain Schrodinger-type partial differential equation and 
proceed to characterise all diffusions which can be interpreted as having the law of B^ conditioned 
on the law of A^ . In the case when oil , . . . , ad is a simple system, these diffusions are closely related 
to the (generalised) quantum Toda lattice. The Schrodinger operator is 

i 

where 9 £ R d , and the corresponding diffusion has generator given by 

-A + VlogVV, 

where fc M is a particular eigenfunction of H known as a class one Whittaker function. We refer to 
these processes as Whittaker processes. In the case n = d — 1 and a\(x) — x, k^{x) — K^(e~ x ) 
and the Whittaker process has generator given by ([lj. In the general 'type A' case, we anticipate 
that the Whittaker process introduced in this paper will have applications, via a multi-dimensional 
generalisation of the theorem of Matsumoto and Yor, to the directed polymer model introduced in 
the paper [22]. This connection is the subject of ongoing research and will be presented in more 
detail in a future work. 

The outline of the paper is as follows. In section 2 we work in a general setting and establish a 
Schrodinger type partial differential equation satisfied by the characteristic function of exponential 
functionals of a multidimensional Brownian motion. We also study a family of martingales related 
to the conditional laws of exponential functionals that will later appear. In section 3, we identify 
a family of diffusions which can be interpreted as having the law of the Brownian motion with 
drift conditioned on the law of its exponential functionals. In section 4, we restrict our attention 
to the case where the collection of vectors used to define the exponential functionals is a simple 
system, and give an overview of relevant facts about class one Whittaker functions. In section 5, 
we study properties of the conditioned processes in this setting. In the final section, we present 
some explicit results for the 'type A^ case. 

2. Laws of exponential functionals and associated partial differential equations 

In this section, we work in a general setting and establish a Schrodinger type partial differential 
equation satisfied by the characteristic function of exponential functionals of a multidimensional 
Brownian motion. We also study a family of martingales related to the conditional laws of expo- 
nential functionals that will later appear. 

Let ai, . . . , ad be a collection of distinct, non-zero vectors in R. n such that fi = {x E R™ : 
at(x) > V?} is non-empty. Let B^) be a standard Brownian motion in W 1 with drift fj, £ f2. For 
< t < oo, set 

A\= f e' 2a ^Us i = l,...,d. 
Jo 

Here, ai((3) = (ai,(3) where (•, •) denotes the usual inner product on W 1 . 

2.1. Partial differential equation for the characteristic function. The process (^B^, A t 
is a diffusion with generator 
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We first check that this operator is hypoclliptic. 
Proposition 2.1. The operator 



2aj(x) a 

da,. 



is hypoelliptic on M. n+d and therefore, for t > the random variable (b^, Atj admits a smooth 
density with respect to the Lebesgue measure. 

PROOF. We use Hormander's theorem. Since the a^'s are pairwise different and non zero, 
there exists d£E" such that 

i + 3 => oti{v) ^ OLj(v). 

Consider now the vector field 

d_ 
dxi 



i=l 

and let us denote 

d 



e -2c«i(x) 



dai 



The Lie bracket between V and T is given by 

d_ 
da,' 



C V T = [V, T] = —2 a,(v)e- 2a ^ - 



i=i 

Similarly, by iterating this bracket k times, we get 

d 



C k v T = (-l) k 2 k a,(v) k e- 2a ' {x) 



d_ 



i=i 

Since the on 's are pairwise different and non zero, we deduce from the Van der Monde determinant 
that at every xel" the family 

{C^T, l<k<d} . 

is a basis of R d . It implies that the Lie bracket generating condition of Hormander is sastisfied so 
that the operator \ A x + (fi, V x ) + Yh=i e~ 2a * (:r) is hypoelliptic. □ 

Let now 6 e R d and, for x e R™, define 

9 l(t,x)=E(e-^^- 2 " i( * )A l), t>0, 

and 

j f l (x)=E{e-^ele-^ A ^y 
Proposition 2.2. (1) The semigroup generated by the Schrddinger operator 



iA + ( M ,V)-]>> 2 e - 2 ^ 

i=l 

admits a heat kernel qt{t,x,y) and we have 

9%x)= [ q%x,y)dy. 
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(2) The function j is the unique bounded function that satisfies the partial differential equa- 
tion 

d 



1*3* {x) + (M, Vj») = £ 6 i e ~ 2ai(X) 3» 



\i=l 



2 

and the limit condition 

lim j e (x) = l. 

x — too.xGil 

Proof. 

(1) It is a straightforward consequence of the Feynman-Kac formula that g^(i,a;,?/) exists 
and is given by 

q l(t,x,y)=E( K e-^ e > e A t \ B y>=y-x)— 

Integrating this with respect to y, we obtain 

gi(t,x)= [ q 6 Jt,x,y)dy. 



(27Ti)t 6 



(2) It is again a straightforward consequence of the Feynman-Kac formula that solves the 
partial differential equation, and the limit condition is easily checked. Let us now prove 
uniqueness. We have to show that if <j> is a bounded solution of the equation that satisfies 

lim 4>{x) = 0, 

x — >co,x£fl 

then = 0. For that, let us observe that under the above conditions, for x £ R™, the 
process 



( f>(B^+x)exp (-J26fe-^Al\ 



is a bounded martingale that goes to when t — > +oo. It follows that this martingale is 
identically zero almost surely, which implies = 0. 

□ 

For later reference, we rephrase the second part of the previous proposition as follows: 
Corollary 2.3. The function ht(x) = e^^jtix) is the unique solution to 



(2) \*h%{x) ~ £ flfe-*"M^(z) = \M 2 h^ 



i=l 



such that e ^ x ^h e ^(x) is bounded and 

lim e-^h 9 Jx) = 1. 

X — >00,3;£Q 

Example 2.1. The following example has been widely studied (see, for example, [8j 118] and 
references therein). Suppose n — d — 1. Q\ — 1/2 and a\(x) = x. Then 

A 00 = [°° e-^+rtdt, m>0 
Jo 

where {B t , t > 0) is a standard one- dimensional Brownian motion, and 

j^^Efexpf-^e-^A^} x£R. 
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d 2 



In this case, h^{x) — e^ x j 6 ^(x) solves the equation 

dx 2 * )K = »K 

This equation is easily solved by means of Bessel functions. By taking into account the boundary 
condition when x — > +00, we recover the formula [18, Theorem 6.2]: 

(3) J» = ^e-^K^e-'), 

where Kn is the Macdonald function |17j : 

1 rx\f f +co e-*-t§ 



(4) K >v=m 



-dt. 



The formula can also be derived using the fact jS] that Aryo has the same law as 1/2^^, where 
7^ is a gamma distributed random variable with parameter fj,. 



Example 2.2. The following example has also been studied in the literature |lll 113) . Suppose 

x_ 
2 ■ 



n = 1, d = 2, 9\ = 6*2 = 1/2, a\{x) — x and a^{x) — §. Then 



coo 

Al= I e^+^dt, At= I e-^+^dt, n>0 



where (B t , t > 0) is a standard one- dimensional Brownian motion, and 
jlix) = E fexp (-le- 2x Al - U-*A a J)) x e 



In this case, hAx) = e^ x jn{%) solves the equation 

d 2 -x -2x\ 1,6 2 u 6 

This is Schrodinger's equation with the so-called Morse potential. It is solved by means of Whittaker 
functions and by taking into account the boundary condition when x — > +00, we get 

j 9 Jx) = 2^ T ^^ e(-^>W 1 J2e~ x ) 

where Wu,^. is the Whittaker function (see [17, pp. 279],): 

^ r + °° „-u»-k-i ( - t^ h h 



WkA*) = wTT TV / e-H^-z 1 + - dt. 

2.2. Conditional densities. We prove now that the random variable Aoo has a smooth den- 
sity with respect to the Lebesgue measure of M. d and morever give an expression of the conditional 
densities only in terms of this density. 

Proposition 2.4. The random variable A^ has a smooth density p with respect to the Lebesgue 
measure of M. d and for t > 

P(4» E dy\F t ) = e^^^p^^Xy,- W)*-*^)^ 
where T is the natural filtration of B^ . 
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Proof. If we denote by <j> the characteristic function of A^: 

0(A) - E ( e -< A ^)) , Ai,...,Ad>0, 

then, 

E (e- £?=i X * A ~ \T t ) = e" Ef=i AiA *E (V Ef =1 |jr t ) 

= e-^)0(e- 2Ql ^" ) )A 1 ,... ,e- 2 ^ s ^A d ) . 

Therefore, the process e~^ A ' A *^ ^e _2ai ^- B *'*^Ai, • • • , e -2 "^- 8 *^) X^j is a martingale. This implies 
that the function e~( A ' a V (e _2ai ( x )Ai, • • • , e~ 2ad ^ \d) is harmonic for the operator ^A x +(fi, V x )+ 
Y?,i=i e_2ai ^gfr- This operator being hypoclliptic, this implies that has a smooth density 
with respect to the Lebesgue measure of R d . The result about the conditional densities stems from 
the injectivity of the Laplace transform. □ 

In particular, we deduce from the previous proposition that if for y e R+, we denote 

q(x,a,y) = e 2 ^ a ^ p (e 2a ^\ Vl - a,),- ■ ■ , e 2a ^(y d - a d j) , < a, < y u x eR d , 

then the process q ,A t , y^j l(o, ! /i)x---x(o,j/ n )( J 4t) is a martingale. It implies that for any y e 
q(a, x, y) satisfies the following partial differential equation: 

2A^+(^,v x g)+]r e - 2 -^ = o. 

1=1 

We can also observe that it implies that p is solution of the partial differential equation: 
Y^(ai,a 3 )+^a t (n) p. 

yi,j=l i=l J 

Example 2.3. Suppose n = d = 1, d\ = 1/2 and ai(x) = x. Then Aoo is distributed as 1/2^^, 
where 7^ is a gamma law with parameter [i, that is 

i_ 



and we have 



. — 2x 



ff( *' °' V) = 2^) \y-a)^ lR >» {y a) 



Example 2.4. Suppose n = 1, d = 2, ot\{x) = x and o^fx) = §. T/ien ? as seen before, 

^ - / e-^'+^dt, ^ = / e-^+^dt, n > 0, 
Jo Jo 
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and for Ai, A2 > 0, 



,r(/'- ! ' 



2 1 2Ai 

r ( 2 ^) -2xr^ 

-t^+2AT 2(2Al+<) Al ^ I (ft. 



r(2Ai) 7 



22 1 



using in the previous integral the change of variable t = e x 13 1 _ 1 , we deduce the following nice 
formula 

^2^+1 Aicotanh^ 

2 (<-~ A > A - 1 .1-; : e ,/« ) = _L_ - 2 ^ +i ds, s >o. 



L-WAlo \A 2 0O = s)¥(A 2 0O €ds) = g—^- 



T/izs conditional Laplace transform can be inverted (see for instance but, unlike the one- 
dimensional case, it does not seem, to lead to a nice formula for p: 

2 2 ' 1 £^ (-iy'2 j T(j + 2u+l + k) 1 g+^+j+M+i)) 2 



r(2 M )V2^.^ j\ fe!r(i + 2/x+i) 

"i + y 2 (fc + j + |) 

where D v is the parabolic cylinder function such that 

r+00 -et 



x D 



j+2/i+2 



/■+00 -et 2 / \ 



that is 

, > \/2 xi /"+ 00 „ _*i / 7TZ/\ , 

A/(z) = -^=e 4 / re 2 CO s I xt — I dt, v > -1 



^ io v 2 y 

3. Brownian motion conditioned on its exponential functionals 

In this section, we study the Doob transforms of the process (^B^\ A t ^j associated with the 
conditioning of A^. We first start with the bridges which are the extremal points. 

Lemma 3.1 (Equation of the bridges). Let y S R+. The law of the process {Bf + /ii)t>o 
conditioned by 

Aoo = y 

solves the following stochastic differential equation: 

dX t = U+(V x \nq) (*t,J e^^da,^ dt + d/3 t 

where, (Pt)t>o is a standard Brownian motion 

Proof. This equation follows directly from Proposition |2.4| and Girsanov's theorem. □ 



Example 3.1. The following example is considered [20) . Suppose n = d = 1, 6\ = 1/2 and 
ct\{x) = x. Then the equation becomes 



e -2X t 
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Let P M be the law of B* 1 and it be the coordinate process on the space of continuous functions 
— > W l . If v is a probability measure on R^., in what follows (see [3j, we call the probability 

+00 r+OO \ 

e- 2a ^ds = »!,•••, / e^^da = y d !/(dy), 



.0 Jo / 

the law of the process (B t + /zt)t>o conditioned by 

Aoo = iat0 1/. 

Proposition 3.1. Let v be a bounded and positive function such that J Rd v(y)p{y)dy — 1. The 
law of the process (B t + //i)t>o conditioned by 

Aoo = law v(x)p(x)dx 

solves the following stochastic differential equation: 

dX t = (a + F v (J e- 2ai{Xs) ds, ■ ■ • , J e - 2ad{x ° ] ds, X^j ^jdt + d(3 t 

where, (flt)t>Q * s a standard Brownian motion and F v : M. d x K™ — * K™ is given by 

F v (a,x) = (V x ki(f> v )(a,x) 

with 

4> v {a,x)= [ p(z)v( ai + e- 2a ^ x) z 1: --- ,a d + e- 2a ^z d )dz. 

Proof. Following [3], we have to write the stochastic differential equation associated with 
the conditioning 

= law p(x)v(x)dx. 

But 

=<t> v (A t ,BW), 

so that we get the expected conditioned stochastic differential equation by Girsanov theorem. 
In the previous proposition, the drift F v (a,x) depends only on x if, and only if, 



□ 



v{x) = e -jm~' 

for some 9 £ M. d . Therefore: 

Corollary 3.2. For 9 £ M. d , the law of the process (B t + p,t)t>o conditioned by 

e~ 



_law 



-p(x)dx 



m 

is the law of a Markov process. Moreover, in that case, it solves in law the following stochastic 
differential equation 

(5) dX t = V In h 8 (X t )dt + d[3 t . 
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We now show that the pathwise uniqueness property holds for the stochastic differential equa- 
tion In what follows, we denote 



£*=Vln/#V + ~A. 



Let us observe that for the generator we have a useful intertwining with the Schrodinger 
operator |A — Yli=i 6fe~ 2oii ^ — ^\\n\\ 2 that will be used several times in the sequel. 



Proposition 3.3. 



Proof. If / is a smooth function then we have 



Since 



2 p 



the result readily follows. □ 
We can now deduce: 

Theorem 3.1. Let 8 £ M. d . If (Pt)t>o is a Brownian motion, then for xq £ M. n , there exists a 
unique process (X^°) t >o adapted to the filtration of (/?t)t>o such that: 

(6) X? = x + f V In hl(Xf°)ds + t , t> 0. 

Jo 

Moreover, in law, the process (X^ a ) t >o is equal to (B t + /it + Xo)t>a conditioned by: 

( / e -2MB t +nt+xo) dt ) = iaw e e ge m, ,e ytf 

\Jo ) l<i< d #(«o) 

Proof. Let xq £ W l . Since the function Vlnft^ is locally Lipschitz, up to an explosion time 
e we have a unique solution X^° for the equation ([6]) . Our goal is now to show that almost surely 
e = +oo. For that, we construct a suitable Lyapunov function for the generator C^. 

Let 

cosh 20,2;) 

U{X) = 

It is easily seen that when || x ||— > +oo, U(x) — > +oo. Moreover, from the intertwining, 



^A -J20?e- 2ai(x) - lM 2 ^J cosh2( M ,a;) < ^|| M || 2 cosh 2( M , x). 



Therefore 

^C/<2 ll/i||2t/ 
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It implies that the process (e - ^"^" t/xe U(X^ e )) t >o is a positive supermartingale. Since U(x) — » 
+00 when || x ||— > +00, we deduce that almost surely e = +00. 

Consequently, there is a unique solution (X x °) t >o for the equation (JsJ) . The second part of the 
theorem is a direct consequence of Corollary |3.2| and uniqueness in law for the equation (|6| . □ 



Example 3.2. Suppose n = d= 1 and Q\ = 1/2 and ui(x) = x. Then 

[i) ^ ^ + e K ^ e - X) ) dx + 2dx 2 

Let us denote p^' (x,y) the heat kernel of C e ^. From the intertwining, we have 

Pt' 9 (x,y) = ^ ' e '_ x , qt' 9 (x,y), 
K^{e x ) 

where q^' e (x, y) is the heat kernel of | — e~ 2x — jJ 2 ^ ■ This kernel can be explicitly computed 
(see pQ or [HI Remark A.l]): 

with 

6(r,t) = —=e™ / e"e- rcoshi; sinh£sin^d£. 

We deduce from that 

p?> e (-ac,y) = 2e-^e(e-y,t)K fl (e-y), 
so that —00 is an entrance point for the diffusion with generator C^. 

The resolvent kernel of (— C e fl + say G^ L,d {x, y, — ^-) is also easily computed: 

G'"' e (x,y,-^-) = 2 ^/^ V \ l r^-^{e- y )K r^-^fe" 31 ), x < y. 
y ' w ' 2 KJe~ x ) V a +c 1 ' V° V ~ 

And we can observe that 

y) = 2K ^~ y ) 1 ^^ 

Example 3.3. Suppose n = 1, d = 2, Q\ = 6 2 = 1/2, a\{x) = x and 0:2(2;) = |. In that case 



and 



re-f 1 2r -, W -^- x A d id 2 
M I 2 W- 1/2!fi (2e-*) ) dx 2dx 2 



u WiJ2e-y) . 
W_x J2e x ) 



where q^ ,S (x, y) is the heat kernel of | (yj^ri — e x — e 2x — . We have (see pQ or |19l p342]j: 



>• + If 



q?' 6 (x, y) = e-^ e -?-(e-'+ e -*)cota„h^ 2 



e-'^r t \ dt, 



sinh £ ' £ / sinh £ 
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G» e (x,y,-—) = — > \ 5 ' M W_, / ^-,(2e-*)M_ 1 -_ 5 - ? (2 e -») ) 



T7ie resolvent kernel of (—C^ + ^-) , is /or x < y; 

r(i + v /^T^) ^_x iM ( 2e -«) 

r(l + 2v/a 2 + /i 2 ) 
and we get: 

a 2 r (1 + t/c^TJA 

2 r (l + 2y/a 2 + ^ 2 J 2 2 ^V«+M 

so £/ia£ — oo is also an entrance point for the diffusion with generator C^. 

Motivated by the two previous examples, the question of existence of entrance laws for the 
diffusion with generator £^ is natural. 
As a general result, we can prove: 

Proposition 3.4. Assume n = 1, ai, • • • , > and 6 e R d - {0}, then — oo is an entrance 
point for the diffusion with generator C ^. 

PROOF. Without loss of generality, we can assume that 9\ > 0. Let us recall h 6 ^ solves the 
Schrodinger equation 

and that k^(x) = K IM {6\e~ aiX ) solves the equation: 

\{1$Y = (fie-***' ** 

Let W{x) = hfc{x){hl)'{x) - (k^Y(x)(hl)(x). Since 

W\x) = k^{x){hl)"{x) - (k^)"(x)(hl)(x) > 0, 
we deduce that W is increasing. Moreover, it is easily seen that lim^^-oo W{x) = 0. Therefore 

W > 0. Hence > -a&e-™ k^'-Z] ■ 

Now, from the comparison principle for stochastic differential equations, we deduce that if, for 
x G R, we denote {Xf ) t >o and (Y t x ) t >o the solutions of the stochastic differential equations, 

t 



o K 



X? = x+ I ^-{X*)ds + f3 t 



where (Pt)t>o is a standard Brownian motion, then we have almost surely 

X? > Y t *. 

Since — oo is an entrance point for the diffusion (Y t x ) t >o, X £R, we deduce that — oo is an entrance 
point for the diffusion with generator £^ . 

□ 
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We conjecture the existence of entrance laws for n > 1, but let us observe that, in general, we 
do not have unicity. Indeed, let us consider the following example 

In that case, by using one dimensional results, we compute: 

pi— ot(p) 

» ( > r(a( M )) A «(^)l e h 

where a*(x) — X2 ^ 1 ■ The heat kernel of is also explicitly given by 

+°° / £ e -2a{x) + e -2a(j/) \ /g-a(x)-Qi(i/) \ ^ 

exp l"2 2^ J (r"c— '*J T' 

And we deduce that when a(a;) — > — oo with a*(cc) -^fegl, 

\J1~Kt 

Therefore, in that case we get an infinite set of entrance laws when a{x) — > — oo. 

4. Whittaker functions 

In this section we restrict our attention to the case when II = {a±, . . . , ad} is a simple system. 
In other words: 

(1) The vectors ot\, . . . , ad are linearly independent; 

(2) the group W generated by reflections through the hyperplanes 

H a = {x G R™ : a(x) = 0}, a £ II 

is finite; 

(3) {x G M. n : a(x) > 0, Va € II} is a fundamental domain for the action of W on W 1 ; 

(4) 2(a, P)/(a, a) G Z for all a, (3 G II. 
In this setting, the Schrodinger operator 

1 d 

is completely integrable and related to the so-called quantum Toda lattice (see [24]). The function 
h e ^ considered in the previous section can be expressed in terms of a particular eigenfunction of H, 
known as class one Whittaker function. 

4.1. Class one Whittaker functions. Class one Whittaker functions associated with semisim- 
ple Lie groups were introduced by Kostant |16j and Jacquet |14j , and have been studied extensively 
in the literature. They are closely related to Whittaker models of class one principal series rep- 
resentations and play an important role in the study of automorphic forms associated with Lie 
groups [6 . In the integrable systems literature, they arise as eigenfunctions of the (generalised) 
quantum Toda lattice |16l 124] . For completeness we will describe briefly the abstract definition 
of class one Whittaker functions, following |12j . Let G be a connected, noncompact, semisimplc 
Lie group with finite centre. Let go be the Lie algebra of G with complexification g. Denote by 
B(-, •) the Killing form on g. Let K be a maximal compact subgroup of G with Lie algebra to 
and denote the complexification of i by t. Let po be the orthogonal complement of to m So with 
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respect to the Killing form. Let 9 be the corresponding Cartan involution. Let do be a maximal 
abelian subspace in p and denote its complexification by a. Denote by S the set of all non-zero 
roots of go relative to do- For denote by m(a) the dimension of the root space 

go = {X G flo : ad(H)X = a(H)X for all H G a }. 

Let S + be a positive system of roots in E and let IT = {a\, . . . ,a.d} be the corresponding set of 
simple roots. Let no = X^es 0o an d N = exp(no). Then G = NAK is an Iwasawa decomposition 
of G. Let -ip be a non-degenerate (unitary) character of N. Let rj be the unique Lie algebra 
homomorphism of no into K such that ip(n) — exp(ir/(X)) for n = cxp(X) G N. For each a G E + , 
let X Qii (1 < i < m(a)) be a basis of satisfying B(X a ^ 7 8X a j) = —Sij (1 < i,j < m(a)). 
Denote by rj a the restriction of T] to Qq and set \r] a \ 2 = J2i<i< m (u) v(X a ,i) 2 - Denote by U(g) and 
U(a) the universal enveloping algebras of g and a, respectively. Let 7 denote the Harish-Chandra 
homomorphism from f7(g) { , the centraliser of t in t/(g), into J7(a). For v G a* and z G C7" (0)* , 
define — 7( z )( iy )- The space of Whittaker functions on G associated with v G a*, denoted 

C^(G/K, Xv), is the space of smooth functions on G which satisfy: 

(1) /(n#fc) = ip{n)f(g) for n e N, g e G and fceif, and 

(2) z/ = X ^)/for2£[/(j) ( . 

Set /? = ^ X Q es + m ( a ) a - For g G G, define ^(g) = h(g)" +p where g — n(g)h(g)k(g) is the 
Iwasawa decomposition of g. Let s be the longest element in W. The cZass one Whittaker 
function associated with v € o* is defined by 

(8) W v {g)= f Uisong^-H^dn, geG. 

J N 

The convergence of this integral was established by Jacquet |14j . For v G a* and a G E, write 
i^a = (a,v)/(a,a). Let 

D = {V G a* : 5R(^ Q ) > 0, for all a G 
We record the following lemma for later reference. 

Lemma 4.1. Let u £ D. Then h~ SaV ~ p W v (h) is uniformly bounded for h G A. 
Proof. Gindikin and Karpelevich |10j proved that the integral 

c{u) = I l u (son)dn 
Jn 

is absolutely convergent. From ^ we can write 

W v (h) = h s ° u+p [ l v (a n)ili~ x (hnh,- l )dn h G A. 
Jn 

Since i/j is unitary, it follows that h~ s ° v ~ p W v {h) is bounded, as required. □ 
Remark 4.1. In the above, c(v) is the Harish-Chandra c-f unction. 

4.2. Fundamental Whittaker functions. Since W u (nhk) — ip(n)W v (h), all of the im- 
portant information about W v is contained in its restriction to A. This leads to a more concrete 
description which can be presented entirely in the context of the root system E. Readers not famil- 
iar with root systems may find it helpful to think of the 'type A' case, for example if G = SL(n, E). 
In this case, we can identify do (and its dual) with 

M' l = {AeR": Ai + • • • + A„ = 0}, 

and take E = {e^ — e,-, i =/= j}, E + = {e^ — e^, 1 < i < j ' < n} and LT = {e^ — e^+i, 2 < i < n}, 
where {e l7 . . . , e„} is the standard basis for M. n . In general, the root system E is crystallographic, 
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that is, the numbers 2(a,f3)/(a,a), a, (3 G H are all integers, and the Z-span of II is a regular 
lattice in oj. Since the Killing form is positive definite on cig, it induces an inner product (•,•) 
on aj, which extends to a nondegenerate bilinear form on o*. The following construction is due 
to Hashizume [12]. Consider the lattice L = 2Z + (LI), and set 'a* = {v G a* : (A, A) + 2(A, v) ^ 
0, VA G L\{0}}. For each v G 'a*, define a set of real numbers {c\{y) : A G L} recursively as follows. 
Set cq(v) = 1 and 

(9) ((A,A)+2(A,^))c A (^) = 2^|7 ?a | 2 CA _ 2Q (^) AgL, 

a 

with the convention that c\{y) = if A ^ L. In |12] it is shown that the series 

AGL 

converges absolutely and uniformly for x G a and ^ G 'a*. Define U to be the set of v G 'a* such 
that: 

(1) u a ^0 for all a G E; 

(2) G 'a* for all s G W; 

(3) sv — tv ^_ X) Q Gn ^ a ^ or an y P a ' r s i'^ such that s 7^ t. 

For s £ W denote by l(s) the length of s. For v £ U, define M(s,v) (s G W), recursively as 
follows. For s — s a (a G II), 

M(s a ,v) = (|77 Q |/2 x /2(a,a))" l/ "e Q (^)e a (-//)- 1 , 

where 

e Q (^)- 1 = r((i/ a + m(a)/2 + l)/2)r((i/ Q + m(a)/2 + m(2a))/2). 
If s G W and a G II such that l(s a s) = l(s) + 1, then 

M{s a s,v) = M(s,v)M(s a , sv). 

Let be the set of a G £+ such that a/2 is not a root. The Harish-Chandra c-function is given 
by 

C{v) = J| d a f a (v), 

where 

= r(i/ )r((i/ + m(a)/2)/2) 

1 J r( I / a + m(a)/2)r((^ a! + m(a)/2 + m(2Q ! ))/2)' 

and 

d Q = 2^"^ a ^" l ' 2a ^/ 2 (7r/(a, a ))( m ( Q ) _m ( 2Q: ))/ 2 . 

Now define, for v G U, 

(10) ¥„(a;) = ^ M(s s,v)c{s sv)® s „(x). 

sew 

Observe that \E'„ satisfies the functional equations 

(11) Vvix) = M(s, u)^ sv {x) s G W. 

Although the above construction places a restriction on v, it is known that, for each x G a, 
^ v {x) can be extended to an entire function of v G a*. In |12j it is shown that, for x G o, 
W v {e~ x ) = e-P^H^x), so that 

W v {g) = i>(n(g))h(gyV„(logh(g)). 
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The functions V v denned by 

V v (g) = ^(n(g))h(gY^(logh(g)), 

are called fundamental Whittaker functions. In [12] it is also shown that, for each v G U, {V sv , s G 
W} form a basis for C$>(G/K, X„). 

4.3. The quantum Toda lattice. As observed by Kostant [16] , Whittaker functions are 
eigenfunctions for the (generalised) quantum Toda lattice. Denote by A the Laplacian on Oo 
corresponding to the Killing form. For v € a*, the class one Whittaker function ^f v (on Oo) 
satisfies the partial differential equation 

(12) \Af(x) - M 2 e- 2a{x) f(x) = \{y,v)f{x). 

For v G U, this can be seen directly via the recursion |9]) for the coefficients in the series expansion 
of the fundamental Whittaker functions $„. In 12 Lemma 7.1] it was shown that, for v 6 D. 

Km e s °" {x) ^ u {x) = c{v). 

x — >oo,xeQ 



By lemma 4.1 if v G D, then e s ° v ^ x ^ v {x) is uniformly bounded for x £ do- Recalling Corol- 
lary |2.3] — note that the proof of uniqueness given there is valid for v G D — we deduce the following 
characterisation of \E' !/ . 

Proposition 4.1. For v G D, the class one Whittaker function ^ v (on <Xq) is the unique 
solution to the partial differential equation 

(13) -A/Or) - \v a \ 2 e- 2a{x) f(x) = \(v,v)f(x), 

such that e s ° v ^ x ^ v (x) is bounded and 

lim e sa " {x) ^ v {x) = c{u). 

x — >oo,xeVL 

4.4. Weyl-invariant class one Whittaker functions and an alternating sum formula. 

In this section we present a variation of the formula ( 10 ) which generalises a formula given in [15 
for the case G = SL(n, R) and leads naturally to a normalisation for the class one Whittaker 
functions which is invariant under the Weyl group W . Using this, we also confirm a conjecture 
of Stade [26 that a class one Whittaker function can be expressed as an alternating sum of 
appropriately normalised fundamental Whittaker functions. 
Let 

a ^= II \ (\v a \/^2{a,a)Y U *T{v a ). 



2 



Proposition 4.2. For v G U, 



c{v) l ^ v {x) = a{v) 1 Y a{s su)^ su {x). 
sew 



PROOF. From (10 1 we have 

= M(s s,v)c(s su)$ sv . 
sew 

It therefore suffices to show that, for all s G W, 

M(s,v)c(sv)a(sv) = c(v)a(v)~ l . 
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We prove this by induction on l(s). If s = s a , (a S II), we have 

M(s a ,v) = {\r la \/2 y /2{a,a) S j 2Va e a {v)e a {-u)-\ 
c{s a v) = / Q (-iy)/ a (^) _1 c(^), 
a{s a y)- l = (\r) a \/^2(a,a)) ^rt^rf-^J-'af.)- 1 . 



Using the duplication formula 

(14) T(z)T(z+ 1 2 )^2 1 -^^T(2z), 
we can write 

and so 



e a {v)/f a {v) = it- 1/2 2^- 1+m ^l 2 /T{v a ) : 
M(s a ,u)c(s a p) = (\i la \/^2{a,a)) 2l ' C ' 



Thus, 



M(s a ,v)c(s a p)a(s a iy) 1 — c{v)a{v) , 
and the claim is proved for l(s) — 1. For s € W and a E II with l(s a s) = l(s) + 1, 

M(s a S, v)c{s a Sv)a{s a Su)^ 1 

= M(s, v)M{s a , sv)c(s a sv)a(s a sv)~ l 
= M{s,v)c{sv)a{svy 1 
= c{v)a{v)~ l , 

by the induction hypothesis. □ 
Consider the normalised Whittaker functions 

w v {x) — a(v)c(v)^ ll i' l ,(x) v e a*, x e a; 

m v (x)= Yl ^\r] a \/^2(a,a)j " T(l + i^)" 1 ^^) v E U, x e a. 
By the above proposition and the functional equation 

(is) r(z)r(i-z) = ^^. 



we have: 

Corollary 4.3. For v e U, 



u v {x)=R{v)- 1 ]T(-l)^W(x), 



sew 

where 

TT 2TO _ 

J.J. 7T 

In particular, w v satisfies the functional equation 

w sv {x) =w v (x), s e W. 

This confirms a conjecture of Stade [26] , who obtained this formula for the case S'L(3,K) and 
conjectured that such a formula holds for all SL(n,M). In the case G = SL(n,M.), the functions 
w v are essentially the same as those considered in [15 . Finally we note that, for each x € a, w l/ (x) 
is an entire function of v. 
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4.5. Plancherel theorem. A Plancherel theorem for class one Whittaker functions was ob- 
tained by Semenov-Tian-Shansky [24 . This is discussed further by Kharchev and Lebedev in [15 
for the case G = SL(n, R). 

PROPOSITION 4.4. The integral transform 

(16) f(x) = J f(T)w lT (x)dT 

is an isomporphism from L 2 (il, (2ir)~ d \a(iT) \ 2 dr) to L 2 (ao,dx), with 

(17) / />)^dz=_I_ [ f( T )g^j\a(zT)\ 2 dT. 

Jao \ Z7T ) JQ 

4.6. The type A\ case. Let G — SX(2,R). Then we can identify ao with R, and take 
E = {±1}, IT = {1} and m(l) = 1. Let |??i| 2 = 1/2. Then L = 2Z+. For A = 2n, write c n = c x {v). 
The recursion (I9j) becomes 4(n 2 + vn)c n = c„_i with cq = 1. The solution is given by 

4" n r(i/+l) 
Cn ~ n!L(n + i/ + 1) ' 

and so 

( p — x /e)\2n+v 

*„(») = mi + v) £ J r(n V, + 1) - 2T(1 + ^-J, 

n>0 v ' 

where I v is the modified Bessel function of the first kind. In this case, W ~ 7L 2 acts on R by 
multiplication. By the duplication formula (14 1, we have 

M( ai ,i/ =4""-) C{U ) = 1 \> 



Thus, using the functional equation (151, we obtain 

* v (aj) = M( Sl , v)c{-v)$ u {x) + c(i/)*_„(a:) = 2 ^- K v (e~ x ), 

where 

7T J_„(z) - J„(z) 

= o ; 

2 sin ttv 

is the Macdonald function. Note that a(A) = 2 A_1 L(A) and the normalised Whittaker functions 
are given by m u {x) — l v (e~ x ) and w v {x) = K v (eT x \ In this setting, the integral transform (16) 
is known as the Kontorovich-Lebedev transform 



4.7. The type A 2 case. Let G = SX(3,R). In this case we can identity ao with Rq = {x G 
R 3 , xi + x 2 +x 3 = 0} and take II = {a\ = (ei - e 2 )/\/2, a 2 = (e 2 - e 3 )/\/2}, where {ei, e 2 , e 3 } is 
the standard basis for R 3 . Set m(ot\) = m(a 2 ) = 1, m{2a\) = m(2a 2 ) = and \rj ai \ 2 = |?7a 2 | 2 = 2. 
For v e R(II) and A = 2nct\ + 2ma 2 e L = 2Z + (II), write c„, m — c\(v). Set a = ct\{v) and 
b = a 2 {v). Then the recursion ^ becomes 

(n 2 + m 2 - nm + an + bm)c n , m = c„_i, m + c„ jm _i, 

where co.o = 1 and c rhm = for (n,m) l\. The solution is given by the following formula, due 
to Bump [6]: 

T(a + l)T(b + l)T(a + b+ l)T(n + m + a + b+l) 
Cn ' m ~ n\m\T(n + a+ l)T(m + b+ l)T(n + a + b + l)T(m + a + b + 1) ' 
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In the notation of [U [7] , 



where 



v 1 = (a+l)/3, i^ = (6+l)/3, yi = 2e- a ^ x \ y 2 = 2e~ a ^ 
The following integral representation is due to Vinogradov and Takhtadzhyan 



(18) 



Wu{x) = -(2/l/y2) " 3 



For a = b = 2/3, we have the following simplification 

2 ' 2 

W(5/9,5/9) (2/1, 2/2j 



V3tt 



K a+b (yWl + r)K a+b (y 2 ^/l + l/r)r~ 
(2/12/2) * (2/| + 2/1 ) 1 * 1/3 ((2/! + 2/1 ) ? 



r 



Using the integral representation (18 1, Bump and Huntley [7, derived an asymptotic expansion of 
W(Vi,«a) (j/i ' 2/2) which is valid for large values of 2/1 and 2/2- The leading term in the expansion is 
independent of the parameter v and given by 



(19) 



/ 2 1 a a _i / a a 

Y ^(2/12/2) 3 (2/1 +2/ 2 3 ) 4 exp(^-(2/ 1 3 + 2/ 2 3 )~ 



From this we deduce the following lemma, which we record for later reference. 
Lemma 4.2. Let A l7 A2 > 0. Ifyi,y 2 — * 00 wzZ/i 2/2/2/1 ~~ * <^> fien 
^W^Cv^ + 2Aij/i, + 2A 22/2) 



IF, 



(^1^2) 



(2/1,2/2) 



exp(-Ai</j(<5) - A 2 <^(l/<5)), 



w/iere 

<p(cf) = (l + di)-di(l + di)i +d^(l + d~i)i. 
4.8. Asymptotics for large x. Consider the analytic function on a* x a defined by 

0(i/,aO=M*)~ 1 E(- 1 ) ,We "' (30 ' 
sew 

where ft,(f) = Ilaes v a- Set 

0* = U(D) = {v £ a* : v a > 0, Va £ n} 

and 

= {.t € a : a(x) > 0, Va e II}. 
PROPOSITION 4.5. Lei q = For all x £fl and v £fl*, 

lim(2c) 9 w_ cl/ (a;/c) = <p(y,x). 

clO 

Proof. First note that, since v £ fi, csv £ U for all s £ W and for all c > sufficiently small. 
The claim follows from Corollary 4.3 and the fact (see |12j ) that there exists a constant k such 
that for all s £ W and c > sufficiently small, 



J2 c x {-c SSoV )e- x ^/ c 
\eL\{o} 



< 



E 



(n + d—l)\ fc" 2mb„ 6 na(i)/c 



J {d-l)\n\ (n!) 



h2 



□ 
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5. Whittaker functions and exponential functionals of Brownian motion 

Define k\ = w_ A for A € a*. Throughout this section we will identify <Xq with ao via the Killing 
form. Let B^' be a Brownian motion in ao with covariance given by the Killing form and drift 
/i G il. Then, by Corollary |2.3| and Proposition |4.1[ we have: 

Proposition 5.1. 



Eexp ( - V \ Va \ 2 e- 2 ^ I e~^ B ^dt ) = e-^c(-s of x)-^^(x) 
= -Q 2 {\r ]a \/^2(a^)y a T(ti a r 1 k fl (x). 



e 



In this context, the diffusion considered in §3 has generator given by 

= + Vlogfc^ • V. 

We will refer to this diffusion is a Whittaker process. Note that it is well-defined for all fj, £ f2. Set 

V,(x) = 2j2\Va\ 2 e- 2a(x) + (^n), 

and write V = Vq. It follows from the intertwining 

(20) k^^^A-V^, 
that the heat semigroup associated with is given by 

(21) P? = e-iW 2 *^^, 

where (Qt) is the heat semigroup associated with |(A — V). By proposition 4.4 the semigroup 
(Qt) is characterised by 

Q t k lT = e"2ll T ll 2 *fc ir , ren. 

Let [j, <E and consider the operator A^, defined (on a suitable domain) by 

A p e A = fc M +A, A € ia* t , 

where e\(x) — e x ^ x \ Set K p = k~ 1 A^. In the type j4i and A 2 cases, for each x, k^ + x(x) is a 
non-negative definite function of A and hence is a Markov operator. For the type Ai case, this 
follows from the integral representation Q and, for the type A^ case, it follows from the integral 
representation (18 1. We conjecture that K M is a Markov operator, in general. In the type A\ case, 
it is shown in [18j that K M intertwines the semigroup associated with with the semigroup of a 
Brownian motion with drift /i. This intertwining extends to the general setting: 

Proposition 5.2. On a suitable domain, 

Proof. For each A £ idg, we have 

(A - V^A^ex = (A - V M )k M+ i 

= (V^+x - V^k^+x 
= ((A,A)+2( M ,A))fc M+A 
= ((A,A) + 2(M,A))A M e A 
= A^(A + 2(i • V)e A . 
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Thus, 



(A-igA„ = A iU (A + 2M-V). 



Combining this with (20 1, we are done. 



□ 



5.1. Brownian motion in a Weyl chamber and Duistermaat-Heckman measure. Let 

ji G il and, for c > 0, define k^(x) = k cjl {xjc). By Proposition 4.5 the diffusion with generator 



££= *A + Vlog^-V, 

converges weakly as c J, to a Brownian motion with drift \x conditioned (in the sense of Doob) 
never to exit the Weyl chamber il (see [4] for a definition of this process). In the limiting case 
/i = 0, the generator of the Brownian moton conditioned never to exit £1 is given by 



-A + Vlog/i- V, 
where h(x) — Yl ae ^o a(x). Note also that, as c J, 0, 

A c e A (x) := (2c)"A Q e cX (x/c) = (2c)«k cX (x/c) - <j>{\,x). 

Thus, the intertwining operator A c converges, in a weak sense, to a positive integral operator with 
kernel given by L(x, dt) = m x DH (dt), where m x DH is the Duistermaat-Heckman measure associated 
with the point igO, characterised by 

e x{t) m x DH (dt) = <p(X,x), Xea*. 

This operator is discussed in [4 . The intertwining 

(A + 2Vlog/i- V)£ = £A 

plays a meaningful role in the multi-dimensional generalisations of Pitman's 2M — X theorem 
obtained in |23[ [5j 2] . 



6. Whittaker processes in the type A 2 case 



Consider the type A 2 case, as in {4.7 For x S Kg, a i( x ) — (x 1 — x 2 )/^/2 and a%(x) = 
(x 2 — x 3 )/\/2. The Weyl chamber is = {x G R 3 , : x 1 > x 2 > x 3 }. Let be a Brownian 
motion in Rq with drift fi £ CI. For < t < oo, set 



i = 1,2. 



Let i/ = — so/i = (—A* 3 , — fJ?, — M 1 )- Then, in the notation of {4.7 

' 3 a+1 

■ v \ = , ^2 



1 2 2 

b = 



V2 



V2 



6 + 1 
3 



Note that, for x G Ttp3 



0; 



= \/2(z/ + ^ 2 )qi(x) + v^f 1 0*2(2;) 

2a + 4b . . 4a + 2b , , 
= — 3 — H 3 — a 2 (x) 

= (2i^x + 4z/ 2 - 2)«i(x) + (4v x + 2v 2 - 2)a 2 (x). 
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By Proposition 5.1 and the integral formula (18 1 



E ( exp y~\yl Al oo-\vl A lc 

cy — a — b 

= **Vi y 2 T(a)T(b)T(a + b) W ^ {y ^ V2) 

o2-a-6 poo j 

= f? 4w I l\ / ^ +b - 1 ^ +& (z/iyr+^)y 2 a+fc - 1 ^ +6 (y 2 v / r+v^)^-- 

r(ajr(&jr(a + 6) J r 

Let us observe that this Laplace transform can be inverted. Indeed, by using the fact that 

p+oo 

K a+b (x)= I e~ xcoshu cosh((a + b)u)du 

we obtain 



I ( exp y-^ylAl, - \y\Al, 

-yi — a—b p+oo p+oo p+oo 

(y 1 y 2 ) a+b ~ 1 / / / fT e~ Vl coshu ^ 1 + r P ~y2 coshtiyl+l/r 



T(a)T(b)T(a + b) Jo Jo 

dv 

x cosh ((a + b)u) cosh((a + b)v) — dudv 

r 

and therefore, denoting p the density of (AJ , ASf), we get 

A p + OO p + OO 

7rr(a)r(6)r(a + 6) J J J 

/ \/i + rcoshtA / ^1 + 1/r coshu \ dr , , 

L> a +b \ ) -IJa+b x cosh( (a + b)u) cosh( (a + o) wl — duaw 



6.1. The intertwining operator. Suppose /i = 0. The intertwining operator A = A satisfes 
Ae_j, = Let a = ai(y), b — a^v) and write t = t±a\ + t 2 a 2 . By the integral representation 



(18) 



2 Jo Jo Jo V uv J \ uv^Jr 

( yWl + r, , 1, 2/2V 1 + !/ r / , 1\ 1 dudvdr 

xexp («+ ~) («+ -) , 

\ 2 u 2 v j u v r 

where 2/1 = 2e~ aityX \ y 2 — 2e~" 2(x ) and 8 = 2/2/2/1 ■ It follows, by a straightforward calculation, 
that we can write 

00 />oo 



Ae- U (x) = / A(x,t)e- atl e~ bt2 dt 1 dt 2 , 

J — 00 J —00 

where 

A(s,t) = K (\fyj{l + Sie^+ t2 )(l + <J§e'i)(l + Sle~^)\ 
A plot of A(x, t), with x fixed, is shown in Figure 1. 
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besselk(0, 1/200 ((1/2+exp(-x+y)) (1/2+exp(-y)) (1/2+exp(x))) 1 ' 2 ) 




y 



x 



Figure 1. The intertwining kernel A(x, •) 

6.2. Behaviour at — oo. From the asymptotic expansion of [TJ, for any A we have k fl (x)~ 1 k\(x) 
1 as x — ► — oo (in the sense that a\(x) — ► — oo and 0:2(2;) — > — oo). This suggests that the Whittaker 
process with generator 

£ M = + Vlogfc^ • V 
has a unique entrance law starting from — oo, given by 

rf{dx) = e-^ 2 %(x)9 t (dx), 
where, for each t > 0, the measure 9t is characterised by 

J k tT {x)9 t {dx) = e -ilMI 2 *, te!1. 

Unfortunately we are unable to provide a rigorous proof of this claim because we do not know if 
the functions k iT are convergence-determining. (The fact that they are determining follows from 
the Plancherel theorem, but this is not quite sufficient to establish uniqueness of the entrance law.) 
On the other hand, we observe: 

Proposition 6.1. Let (X^°) t >o be the diffusion with generator started at xq. If a\(xo) — > 
— oo and a2(x ) — > — oo, with ai(xo) — ^(^o) ~~ * K > then 

/ F + OO /-+CX0 \ 

r e ai(«o) J e - 2a ^ x ^ds,e a ^ J e~ 2a - {x ^ds\ 

converges in probability to (tp(e K ),ip(e~ K )), where 

ip(d) = (1 + di) -dt(i + df)s +d3(l + d-§)i. 
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Proof. Let Ai,A2 > 0. We easily compute 

<>+oo 



E exp -A/ 1 * 10 ' 



+ 00 



e -2a 2 (x;o )ds 



= Eexp( -|(yg + 2X 1 y 1 )A 1 00 - \{y\ + 2X m )A^) 
with 2/1 = 2e- Ql ( :Eo ) 



Ee X p(-iyMio-|^ 2 oo) 
,2/2 = 2e — a2 ( a °). But, by Lemma 4.2 if 2/1, 2/2 — > 00 with 2/2/2/1 - * 5 

. _-A lV (5)-A 3¥3 (l/5) 



Eexp (-|(2/f + 2Ai2/i)^ - |(2/ 2 2 + 2A 2 j/ 2 )^ 



Eexp(-i2/2AL-i2/ 2 2 ^) 



then 



□ 
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